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At r act--  Falling  of  mine  through  air,  water,  and 
sediment  is  investigated  experimentally  and  theoretically. 
Two  experiments  were  conducted  to  drop  cylindrical  mine 
with  the  density  ratio  around  1.8  into  shallow  water  (around 
13  m  deep)  in  the  Monterey  Bay  (Exp-1)  and  into  the  Naval 
Postgraduate  School’s  swimming  pool  (Exp-2).  During  the 
experiments,  we  carefully  observe  mine  track  and  burial 
depth  while  simultaneously  taking  gravity  cores  (in  Exp-1). 
After  analyzing  the  gravity  cores,  we  obtain  the  bottom 
sediment  density  and  shear  strength  profiles.  The 
theoretical  work  includes  the  development  of  3D  mine 
impact  burial  prediction  model  (IMPACT35)  which 
contains  three  components:  triple  coordinate  transform, 
hydrodynamics  of  falling  rigid  object  in  a  single  medium 
(air,  water,  or  sediment)  and  in  multiple  media  (air-water 
and  water-sediment  interfaces).  The  model  predicts  the 
rigid  body’s  trajectory  in  the  water  column  and  burial  depth 
and  orientation  in  the  sediment.  The  experimental  data 
(burial  depth,  sediment  density  and  shear  strength)  are  used 
to  evaluate  the  newly  developed  numerical  model.  The  3D 
model  shows  great  improvement  to  the  currently  used  US 
Navy’s  2D  model  (i.e.,  IMPACT28). 

1.  Introduction 

The  conclusion  of  the  cold  war  culminated  with  the  Union 
of  Soviet  Socialist  Republics  (USSR)  effectively  ceasing  to 
exist  under  international  law  on  December  31,  1991.  This 
historical  event  caused  the  U.S.  military  and  specifically  the 
Navy  and  Marine  Corp  Team  to  shift  tactical  emphasis 
from  blue  water,  deep  ocean  doctrine  to  littoral  warfare 
doctrine.  This  shift  predicated  military  responses  dealing 
with  a  wide  range  of  worldwide  regional  crises  requiring 
forward  sea  basing,  and  expeditionary  force  landing 
support. 

The  Navy  Marine  Corp  team  developed  a  doctrine  concept 
white  paper,  “...  From  the  Sea,  1992”,  to  support  joint 

Manuscript  received  January  15,  2004.  This  work  was  supported  by  the 
U.S.  Office  of  Naval  Research,  Naval  Oceanographic  Office,  and  Naval 
Postgraduate  School. 

Peter  C.  Chu  is  with  the  Naval  Postgraduate  School,  Monterey,  CA 
93943  USA  (telephone:  831-677-3688,  e-mail:  chu@  nps.navy.mil). 


warfare  doctrine  concepts  of  forward  presence  and 
engagement  developed  as  National  Defense  Strategy  in 
“Joint  Vision  2010,  1996”,  (Rhodes  and  Holder  1998).  The 
document  provided  guiding  tenets  for  naval  operations  of 
the  21st  century.  A  subsequent  Naval  Department  revision, 
“Forward  ...From  the  Sea,  1994”,  and  its  Marine  Corp 
counterparts  “Operational  Maneuver  from  the  Sea,  1996”, 
and  “Ship  to  Objective  Maneuver,  1997”  all  focus  on  sea 
based  power  projection  into  littoral  regions  and  guiding 
naval  operations  in  those  areas  in  the  new  millennium. 
“Joint  Vision  2020,  2000”  and  “Sea  Strike,  Sea  Shield,  Sea 
Basing,  2002”  are  the  current  National  Defense  Strategy 
and  Naval  Department  concept  papers  providing  guiding 
tenets  for  naval  and  joint  operations  well  into  the  21st 
century.  Both  papers  incorporate  emerging  technology, 
processes,  people  and  organizations  synergized  via  the 
netcentric  warfare  concept  to  provide  total  power  projection 
and  dominance  across  littoral  regions  during  any  crises 
requiring  U.S.  response. 

...the  very  shallow  water  (VSW)  region  is 
a  critical  point  for  our  offensive  forces 
and  can  easily,  quickly  and  cheaply  be 
exploited  by  the  enemy.  The  magnitude 
of  the  current  deficiency  in 
reconnaissance  and  neutralization  in  these 
regions  and  the  impact  on  amphibious 
assault  operations  were  demonstrated 
during  Operation  Desert  Storm.  Maj. 

Gen.  Edward  J.  Hanlon  Jr.  “From 
(Rhodes  and  Holder  1998).” 

Any  military  operation  that  occurs  in  the  littoral  regions 
also  occurs  in  mine  country.  The  increasing  pace  of 
shallow- water  naval  operations  (i.e.  Persian  Gulf,  Adriatic 
Sea,  Yellow  Sea,  and  Gulf  of  Aden)  translates  into  a  high 
probability  of  encountering  mines.  The  required  shift  in 
focus  of  naval  operations  from  the  open  ocean  to  the 
regional  littoral  areas  increases  the  importance  of  mine 
warfare  as  a  navy  core  competency.  The  proliferation  of 
inexpensive,  bottom  type  mines  make  shallow  water  and 
very  shallow  water  MCM  a  critical  and  expensive 
challenge.  In  times  of  conflict  domination  of  coastal 
operating  areas  will  largely  depend  on  the  ability  to  remove 
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or  neutralize  any  emplaced  littoral  mine  threat,  figure  1, 
and  prepare  the  battle  space  for  follow-on  action  in  a  timely 
fashion  Naval  mines  may  be  found  throughout  the  water 
column  and  on  or  within  the  seafloor.  Ask  anyone  to 
describe  a  typical  naval  mine  to  you  and  the  response  will 
be  a  description  of  the  spherical,  hertz-horn  World  War  II 
vintage  drifting  mine  shape  common  in  Hollywood  films. 
But,  it  is  the  buried  naval  mine  that  poses  the  most  severe 
threat  to  naval  assets  since  naval  forces  possess  very 
limited  resources  and  capabilities  for  detecting,  identifying, 
and  neutralizing  them,  and  the  mine  itself  remains  fully 
effective  when  buried,  (Lott,  2001).  An  important  factor  in 
mine  hunting  and  clearance  is  the  amount  of  initial  impact 
and  subsequent  sediment  burial  a  mine  undergoes  with  time 
because  buried  mines  are  substantially  more  difficult  to 
detect  and  classify.  The  amount  of  burial  becomes  a  critical 
parameter  and  crossroad  in  the  naval  MCM  mission 
planning  process  because  the  mine  countermeasures  effort 
transitions  from  mine  hunting  to  mine  sweeping,  (Rennie 
2002.)  Littoral  Mine  Threat.  “From  (Rhodes  1998).” 

Mine  warfare,  perhaps  more  than  any 
other  single  littoral  warfare  mission  area 
is  the  “key”  that  will  unlock  the  “door”  to 
the  littoral  battle  space.  In  the  most 
fundamental  way,  then,  mine  warfare  and 
the  need  for  effective  mine 
countermeasures  must  be  an  “all-hands” 
concern  for  the  Navy  and  the  Marine 
Corps.  (Boorda  1995) 


Fig.l.  Littoral  Zone  Mine  Treat  (from 
Rhodes  1998). 

The  Office  of  Naval  Research  (ONR)  in  1999  created  the 
Mine  Burial  Program  (MBP),  an  applied  (6.2)  research 
program  to  develop  mine  burial  prediction  tools  and 
decision  aids.  The  program  mission  is  to  predict  the 
behavior  of  mines  in  different  environments,  (Bennett 
2000).  The  Impact  Mine  Burial  Prediction  model 
development  falls  under  the  MBP  program. 

The  Impact  Mine  Burial  Prediction  model  developed  in 
1980  was  developed  to  semi-empirically  model  the  mine 


burial  process.  Several  revisions  have  occurred  in  the  last 
two  decades  but  there  have  been  limitations  noted  in  the 
model  performance,  as  well  as  little  scientific  advancement 
in  mine  burial  prediction,  (Dolan  et  al  1999),  (Taber  1999), 
(Smith  2000),  and  (Gilless  2001). 

In  this  study,  a  nonlinear  dynamical  system  (IMPACT35) 
is  established  for  the  movement  of  a  nonuniform  (center  of 
gravity  not  the  same  as  the  center  of  volume)  mine  through 
the  water-sediment  interface.  A  cylinder-drop  experiment 
was  conducted.  The  data  collected  from  the  experiment  can 
be  used  for  model  development  and  verification.  The  long¬ 
term  goal  is  for  model  inclusion  in  a  full  spectrum 
deterministic  model  for  mine  burial  in  a  comprehensive 
navy  tactical  decision  aid. 

2.  Triple  Coordinate  Systems 

Consider  an  axially  symmetric  cylinder  with  the  centers 
of  mass  (COM)  X  [or  called  gravity  center  (GC)  in 
literatures]  and  center  of  volume  (COV)  B  on  the  main  axis 
(Fig.  2a).  Let  (Z,  R,  %)  represent  the  cylinder’s  length, 

radius,  and  the  distance  between  the  two  points  (X,  B).  The 
positive  X  -values  refer  to  nose-down  case,  i.e.,  the  point  X 

is  lower  than  the  point  B.  Three  coordinate  systems  are 
used  to  model  the  falling  cylinder  through  the  air,  water, 
and  sediment  phases:  earth-fixed  coordinate  (E-coordinate), 
main-axis  following  coordinate  (M-coordinate),  and  force 
following  coordinate  (F-coordinate)  systems.  All  the 
systems  are  three-dimensional,  orthogonal,  and  right- 
handed  ([2]  Chu  et  al.  2004). 

The  E-coordinate  is  represented  by  FE(0,  i,  j,  k)  with  the 
origin  ‘O’,  and  three  axes:  x-,  y-  axes  (horizontal)  with  the 
unit  vectors  (i,  j)  and  z-axis  (vertical)  with  the  unit  vector  k 
(upward  positive).  The  position  of  the  cylinder  is 
represented  by  the  COM  position, 

X  =  xi+yj+zk.  (1) 

The  translation  velocity  is  given  by 
dX 

—  =  V,  V  =  (w,  v,  w).  (2) 

dt 

Let  orientation  of  the  cylinder’s  main-axis  (pointing 
downward)  is  given  by  iM,  and  ^/jbe  the  angle  that  the 
cylinder  rotates  around  the  axis  iM.  The  angle  between  iM 
and  k  is  denoted  by  y/2  +  71  /  2  .  Here,  y/2  is  usually  called 
the  elevation  angle.  Projection  of  the  vector  iM  onto  the  (x, 
y)  plane  creates  angle  between  the  projection  and  the 

x-axis  (Fig.  2b).  The  M-coordinate  is  represented  by  FM(X, 
iM,  )m?  kM)  with  the  origin  ‘X’,  unit  vectors  (iM,  kM),  and 
coordinates  (x(m\  y(m\  z(m)).  The  unit  vectors  of  the  M- 
coordinate  system  are  given  by  (Fig.  2b) 

J m  =  k  x  iM ,  kM  =  iM  x  jM  (3) 
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The  M-coordinate  system  is  solely  determined  by 
orientation  of  the  cylinder’s  main-axis  iM. 

The  F-coordinate  is  represented  by  Ff(X,  iF,  jF,  kF)  with 
the  origin  X,  unit  vectors  (iF,  jF,  kF),  and  coordinates  (x(F) , 
z™).  Let  Vw  be  the  fluid  velocity.  The  water- to-cylinder 
velocity  is  represented  by 

Vr  =  Vw-V,  (4) 

which  can  be  decomposed  into  two  parts, 

V  =V|+V2,V|=  (Vr  ■  iF)iF ,  v2  =  V,.  -  (V  •  iF  )iF  (5) 

where  Vi  is  the  component  paralleling  to  the  cylinder’s 
main-axis  (i.e.,  along  iM),  and  V2  is  the  component 
perpendicular  to  the  cylinder’s  main-axial  direction.  The 
unit  vectors  for  the  F-coordinate  are  defined  by  (column 
vectors) 

J*f=V2/|V2|,  kF=iFxjF.  (6) 


3.  Dynamics 

3.1.  Momentum  Balance 

The  translation  velocity  of  the  cylinder  (V)  is  governed 
by  the  momentum  equation  in  the  E-coordinate  system, 


u 

"o" 

V 

=  - 

0 

w 

_g_ 

F.+F/, 

pU 


(7) 


where  g  is  the  gravitational  acceleration;  Fb  is  the  buoyancy 
force;  n  is  the  cylinder  volume;  p  is  the  rigid  body 

density;  p  n  =  m,  is  the  cylinder  mass;  is  the 
hydrodynamic  force  (i.e.,  surface  force  including  drag,  lift, 
impact  forces).  The  drag  and  lift  forces  are  calculated 
using  the  drag  and  lift  laws  with  the  given  water-to-cylinder 
velocity  (Vr).  In  the  F-coordinate,  Vr  is  decomposed  into 
along-cylinder  (Vi)  and  across-cylinder  (V2)  components. 


3.2.  Moment  of  Momentum  Equation 

It  is  convenient  to  write  the  moment  of  momentum 
equation 

d(0 

J - =  M6+Ma,  (8) 

dt 

in  the  M-coordinate  system  with  the  cylinder’s  angular 
velocity  components  (co1,co2,co3)  defined  by  (4).  The 

gravity  force,  passing  the  COM,  doesn’t  induce  the 
moment.  M^  and  M^  are  the  buoyancy  and  hydrodynamic 
force  torques.  In  the  M-coordinate  system,  the  moment  of 
gyration  tensor  for  the  axially  symmetric  cylinder  is  a 
diagonal  matrix 


J  = 


J, 

0 

0 


0 

J, 

0 


0 

0 

^3_ 


(9) 


where  J\,  J2,  and  J3  are  the  moments  of  inertia.  The  gravity 
force,  passing  the  center  of  mass,  doesn’t  induce  the 
moment.  The  buoyancy  force  induces  the  moment  in  the  jM 
direction  if  the  COM  doesn’t  coincide  with  the  COV  (i.e., 
X*  0), 


M6  =|Fikc°s^2jM-  (10) 

Computation  of  buoyancy  and  hydrodynamic  forces  (Fb, 
Fb)  and  torques  (Mb,  Mb)  is  more  complicated  for  a 
cylinder  penetrating  through  air-water  and  water-sediment 
interfaces  than  falling  through  a  single  medium  such  as 
water.  At  the  instance  when  the  cylinder  penetrates  into  an 
interface,  three  situations  may  exist:  one-side  entry  with 
the  lower  surface  partially  contacting  the  interface  (Fig. 
3a),  one-side  entry  with  partial  area  with  the  lower  surface 
completely  inside  the  lower  medium  (Fig.  3b),  and  two-side 
entry  (Fig.  3c).  For  the  two- side  penetration,  the  upper 
(above  the  interface)  and  lower  (below  the  interface)  parts 
of  the  cylinder  are  represented  by  D(1)  and  D(2).  For  the  one- 
side  penetration,  the  upper  (or  lower)  part  of  the  cylinder 
are  considered  as  the  combination  of  D(1)  [or  D(2)]  and  a 
sub-cylinder  [or  C52)].  In  general,  the  upper  and  lower 
parts  of  the  cylinder  are  represented  by  [C^l\  D(1)]  and  [CP\ 
D(2)].  For  the  two-side  penetration,  =  CP=  0. 


4.  Buoyancy  Force 

Let  ( Lc>  Lm )  be  the  lengths  of  [C^,  Z)(1)],  (hi,  h2 )  the 
heights  of  the  two  sides  of  D(1),  ( F^  ,  F^2) )  the  buoyancy 


forces  for  the  upper  and  lower  parts  of  the  cylinder  (Fig.  4). 
The  vertical  cross-section  of  D is  shown  in  Fig.  5.  In  the 
M-coordinate  system,  we  have 


/  (m)\  D2  “1 

s(x  )  =  R  cos 


f 


1-^ 


Ah 


LR 


J» 


\ 

J 


(ii) 


where  Ah  =  h2  -  \ .  The  integration  of  s(x)  along  x(w)  axis 
gives  the  volume  of  D(1), 


n(1)  =  [  m  s(x(w))dx(w)  = - —  P(kx,k2)  , 

Jo  Ah 

where 

I -  1  - 

^2)  =  kx  cos-1  (k-j )  -  VI  -  K\  +  - 11  -  K\  )2 


-1 

-k2  cos 


3 


(12) 
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kx  -  1  —  hx  /  R,  k2  =  1  -  h2  /  R  . 
The  volume  of  is  calculated  by 


Q 


(i) 


kR2Lc. 


(13) 


Fig.  2b.  Three  coordinate  systems. 


(14) 


The  buoyancy  force  is  the  product  of  the  medium  (air, 
water,  or  sediment)  density  and  volume  fl(1)  +  , 


F^(1)  =  p(l\u(l)  +Qvi;)k 


ON 


JO 


Ah 


P(kvk2)  +  7iR  L 


(15) 


The  volume  of  the  lower  part  of  the  cylinder  [C*2),  ZJ(2)]  is 
given  by 

n(2)  +  Q(2)  =  ttR2L  -  n(1)  -  Q(1) .  (16) 


fe) 


Fig.  3.  (a)  One-side  entry  with  the  lower  surface 
partially  contacting  the  interface,  (b)  one-side  entry 
with  partial  area  with  the  lower  surface  completely 
inside  the  lower  medium,  and  (c)  two-side  entry. 


(2) 

The  buoyancy  force  F^  is  calculated  by 

Ff}  =  p(2)(n(2)  +Q(2))k. 

Substitution  of  (15)  and  (16)  into  (17)  leads  to 


f i2)=p(2) 


7iR  (L-L  )-- 


Ah 


-P(kvk2) 


(17) 


k.  (18) 


(a) 


(b) 


J  m 


Fig.  2a.  M-coordinate  with  the  COM  as  the  origin  X 
and  (im,  jm)  as  the  two  axes.  Here,  X  *s  distance 

between  the  COV  (B)  and  COM,  (Z,  R )  are  the 
cylinder’s  length  and  radius. 


Fig.  4.  Equivalent  cylinder  method  for  computing 
buoyancy  force  and  torque  of  the  volume  above  the 
interface. 


5.  Torque  due  to  the  Buoyancy  Force 

The  torque  due  to  the  buoyancy  force  should  be 
calculated  separately  for  the  upper  and  lower  parts  of  the 
cylinder  as  it  penetrates  the  interface.  Here,  we  use  the 

upper  part  as  the  illustration.  Let  be  the  COVi 

location  of  the  upper  part  of  the  cylinder  in  the  M- 
coordinate  system.  Note  that  the  origin  of  the  M-coordinate 
is  the  COM  of  the  cylinder.  The  location  of  COVi  is 
defined  by 


4 


Journal  of  Counter-Ordnance  Technology  (Sixth  International  Symposium  on  Technology  and  Mine  Problem,  NPS,  10-13 
May,  2004) 


5 


(x 


(i) 

B 


1 

n(1)  +  q(1) 


,z(m))ddX) 


The  torque  due  to  the  buoyancy  force  for  the  lower  part  of 
the  cylinder  can  be  calculated  using  the  same  procedure. 

6.  Hydrodynamic  Force  and  Torque  for  Air 
and  Water 


1 


n(1)  +  q0) 


r  (xM,zM)s(xM)dxM  +  (^,0)d 
Jo  2 


(19) 


Substitution  of  (1 1),  (12)  and  (14)  into  (19)  leads  to 


JO 


RAh 


p(KrKi)Lm(\+nAhLjr't) 


L 
Ah 

2\R 


,(i) 


R 
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where 


1 


(2 k\  -  l)cOS  1  K2  -  {2^  -  l)cOS  1  Kl 


+kJ\-kx  —  k2\\\  —  k. 


1 

+  — 
4 


x\ 


.5 

3  L 


v 

/Z  ; 

s  V 

)  - 

) 

r.  r  . 

-1 

1 1  -  kx  +  sm  b 

-  Kx  COS  1  Kx  +  yj 

r - 

1-/ 

'-v 

? 

+  —  (aTj  ^1  -  ATj2  -  ?c2  aJi  -  /c22  +  sin-1  at,  -  sin-1  a:2  j 


//;(a-,a-2)  =  a- 

3 
2 


(22) 


The  distance  between  COM  and  CO V 1  is  calculated  by 

TU)  =  •  (23) 

The  torque  due  to  the  buoyancy  force  for  the  upper  part  is 
given  by  [see  (10] 

(24) 


M 


(i)  _ 


T~i(l )  (1) 

Fb  z  cos^2jM. 


Substitution  of  (18)  and  (23)  into  (24)  leads  to 
P(kvk2) 


M(d  =  pW 


r3l 


Ah 
+xR2L 


L  (1)  -.2  r  (1)  -.2 

VtXfl  1  +  \-ZB  1  C0S^2JM- 

(25) 


Let  (cr(1),cr(2))  be  the  surface  of  the  upper  and  lower 
parts  of  the  cylinder.  During  the  penetration,  the  total 
hydrodynamic  force  and  associated  torque  are  decomposed 
into  upper  and  lower  parts, 

F,  =  +  Ff } ,  M*  =  M<1}  +  ,  (26) 

where  (F^M^)  and  (F^2),M^2) )  are  the  force  and 

torque  on  the  upper  and  lower  surfaces  a ^  and  cr(2)  . 
Let  the  hydrodynamic  force  be  fh  at  a  point  r  (represented 
in  M-coordinate)  on  the  cylinder’s  surface,  the 
hydrodynamic  forces  and  torques  are  calculated  by 

FA(1)  =  J1  fhdc L  Ff )  =  11  f \d(7’ 

cr,  cr? 

\  ,  •  (27) 

Mi1>  =  11  (r  x  f h)d(J’  Mi2)  =  11  (r  x  h)da 
G\  G2 

The  hydrodynamic  force  and  torque  for  air  and  water  are 
calculated  using  the  recursive  model  recently  developed 
([2]  Chu  et  ah,  2004).  However,  the  surface  force  and 
torque  due  to  sediment  should  be  computed  separately. 

7.  Hydrodynamic  Force  and  Torque  for 
Sediment 

Let  v  be  the  velocity  at  point  r  on  the  cylinder  surface, 
v  =  V  +  co  x  r  . 

The  velocity  (v)  and  hydrodynamic  force  due  to  the 
sediment  ( f^sed) )  have  normal  and  tangential  components 

(Fig.  6), 

v  =  v„+vr>  v„  =  (v.n)n, 

h  ~K  +  ’  (2®) 

where  n  is  the  normal  unit  vector  (outward  positive);  vn 
and  \T  are  normal  and  tangential  components  of  v;  f^sed^ 

and  fr( sed ^  are  the  resistant  forces  in  the  normal  and 
tangential  directions.  Note  that  the  unit  tangential  vector  T 
is  in  the  opposite  direction  of  \T  .  Let  B( z)  and  S(z)  be  the 

sediment  bulk  density  and  shear  strength  profiles.  As  the 
cylinder  penetrates  into  the  sediment,  the  impact  (resistant) 
force  exerted  on  the  part  of  the  object’s  surface  moving 
towards  the  sediment  (Fig.  7), 

f„(W)(z)  =  -n{f;  [B(z')-pwW+S(z)}S,  (29) 

Jzws 
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where  pw  is  water  density  and  8  is  a  step  function 


1  vn  >  0 

S  =  \  •  (30) 

[0  vn  <  0 

Since  the  sediment  is  quasi-solid,  the  resistant  force  in  the 
tangent  direction  (f^e^)  may  be  treated  as  the  friction 


between  two  solid  bodies  that  is  proportional  to  the 
resistant  force  in  the  normal  direction, 


Med)  _  r  Med ) 

~  ln  ’ 


(31) 


where  C  is  the  friction  coefficient  (0.01)  between  the 
cylinder  and  sediment.  The  total  surface  force  due  to 
sediment  is  calculated  by 


j^(sed) 

bh 


+  f(sed) 


M  [sed) 


jj  [rx(fnM)+f (Tsed))da, 


(32) 


®sed 


where  o sed  is  the  area  of  the  cylinder’s  surface  contacting 
with  the  sediment. 


8.  Model  Integration 


The  momentum  equation  (7)  and  moment  of  momentum 
equation  (8)  are  integrated  numerically  using  the  triple 
coordinate  transformation  ([2]  Chu  et  al.  2004).  As  the 
cylinder  penetrates  into  the  sediment,  the  resistant  force  due 

to  sediment  f^sed^  reduces  the  cylinder’s  speed  and  turning 
angle.  The  vertical  coordinate  of  COM  [i.e.,  z(t)]  in  the  E- 
coordinate  and  the  elevation  angle  y/2{t)  are  used  to 


identify  the  cylinder’s  final  burial  depth  and  orientation, 


dz  dii/, 

—  =  0,  =  0  . 
dt  dt 


(33) 


Fig.  6.  Momentum  and  angular  momentum  balance  for 
cylinder’s  penetration  through  the  water-sediment 
interface. 


Fig.  7.  The  impact  (resistant)  force  exerted  on  the  part 
of  the  object’s  surface  moving  towards  the  sediment. 

9.  Cylinder  Drop  Experiments 

Two  cylinder  drop  experiments  were  conducted  to 
collect  data  for  the  model  evaluation.  Exp-1  was  designed 
to  collect  data  on  cylinder's  motion  in  the  water  column  for 
various  combinations  of  the  cylinder’s  parameters.  Exp-2 
was  designed  to  collect  synchronized  data  on  sediment 
parameters  (shear  strength  and  density)  and  the  cylinder’s 
burial  depth  and  orientation. 

9.1.  Exp-1 

Exp-1  was  conducted  at  the  NPS  swim  pool  in  June 
2001.  It  consisted  of  dropping  each  of  three  model 
cylinders  (Fig.  8)  into  the  water  where  each  drop  was 
recorded  underwater  from  two  viewpoints.  The  physical 
parameters  of  the  model  cylinders  are  listed  in  Table  1.  Fig. 
9  depicts  the  overall  setup.  The  controlled  parameters  for 
each  drop  were:  L/R  ratio,  X  -value,  initial  velocity  (Vin), 
and  drop  angle.  The  E-coordinate  system  is  chosen  with  the 
origin  at  the  comer  of  the  swimming  pool  with  the  two 
sides  as  x-  and  y-axes  and  the  vertical  z-axis.  The  initial 
injection  of  cylinders  was  in  the  (y,  z)  plane  (Fig.  10). 

Initial  velocity  (Vin)  was  calculated  by  using  the  voltage 
return  of  an  infrared  photo  detector  located  at  the  base  of 
the  cylinder  injector.  The  infrared  sensor  produced  a  square 
wave  pulse  when  no  light  was  detected  due  to  blockage 
caused  by  the  cylinder's  passage.  The  length  of  the  square 
wave  pulse  was  converted  into  time  by  using  a  universal 
counter.  Dividing  the  cylinder's  length  by  the  universal 
counter's  time  yielded  Vin.  The  cylinders  were  dropped 
from  several  positions  within  the  injector  mechanism  in 
order  to  produce  a  range  of  Vin.  The  method  used  to 
determine  Vin  required  that  the  infrared  light  sensor  be 
located  above  the  water's  surface.  This  distance  was  held 
fixed  throughout  the  experiment  at  10  cm. 
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Fig.  8.  Internal  components  of  the  model  cylinder. 


Fig.  9.  Exp-1  equipments. 


The  drop  angle  (initial  value  of  y/^)  was  controlled 

using  the  drop  angle  device.  Five  screw  positions  marked 
the  15°,  30°,  45°,  60°,  and  75°.  The  drop  angles  were 
determined  from  the  lay  of  the  pool  walkway,  which  was 
assumed  to  be  parallel  to  the  water's  surface.  A  range  of 
drop  angles  was  chosen  to  represent  the  various  entry 
angles  that  air  and  surface  laid  mines  exhibit  in  naval 
operation.  This  range  produced  velocities  whose  horizontal 
and  vertical  components  varied  in  magnitude.  This  allowed 
for  comparison  of  cylinder  trajectory  sensitivity  with  the 
varying  velocity  components. 


For  each  drop  the  cylinder  was  set  to  a  X  -value.  For 
positive  x  -value,  the  cylinders  were  placed  into  the 

injector  so  that  the  COM  was  located  below  the  geometric 
center.  For  negative  X  -value,  the  COM  was  located  above 

the  geometric  center  to  release.  A  series  of  drops  were  then 
conducted  in  order  of  decreasing  mine  length  for  each 
angle.  Table  2  indicates  number  of  drops  conducted  for 
different  drop  angles  and  X  -value  for  L/R  =  15/2.  Number 

of  drops  for  other  L/R  ratios  (12/2,  9/2)  is  comparable  to 
that  for  L/R  ratio  of  15/2.  All  together  there  were  712 
drops.  Each  video  camera  had  a  film  time  of  approximately 
one  hour.  At  the  end  of  the  day,  the  tapes  were  replayed  in 
order  to  determine  clarity  and  optimum  camera  position. 

Upon  completion  of  the  drop  phase,  the  video  from  each 
camera  was  converted  to  digital  format.  The  digital  video 
for  each  view  was  then  analyzed  frame  by  frame  (30  Hz)  in 
order  to  determine  the  mine's  position  in  x-z  and  y-z  planes. 
The  cylinder’s  top  and  bottom  positions  were  input  into  a 
MATLAB  generated  grid,  similar  to  the  ones  within  the 
pool.  The  first  point  to  impact  the  water  was  always  plotted 
first.  This  facilitated  tracking  of  the  initial  entry  point 
throughout  the  water  column.  The  cameras  were  not  time 
synchronized;  thus,  the  first  recorded  position  corresponded 
to  when  the  full  length  of  the  mine  was  in  view. 

9.2.  Exp-2 

Exp-2  was  conducted  on  the  R/V  John  Martin  on  May 
23,  2000.  The  barrel  with  density  ratio  of  1.8  was  released 
while  touching  the  surface.  This  would  be  to  eliminate  any 
chance  of  inertial  effects  caused  by  uneven  introduction 
into  the  air-sea  interface.  This  also  set  the  initial  velocity 
parameter  in  the  code  to  zero.  The  barrel  was  to  be  released 
17  times.  The  diver  would  snap  the  quick-release  shackle 
on  the  barrel  and  then  dive  down  to  conduct  measurements. 
The  average  depth  of  the  water  was  13  meters.  Since  it  was 
uncertain  the  path  the  barrel  would  follow,  both  the 
releasing  diver  and  a  second  safety  diver  would  stay  on  the 
surface  until  after  the  barrel  had  dropped.  Once  reaching 
the  bottom,  one  diver  would  take  penetration  measurements 
using  a  meter  stick  marked  at  millimeter  increments  while 
the  other  would  take  a  gravity  core.  After  17  drops,  the 
divers  began  to  run  out  of  air  and  results  were  not  varying 
greatly  so  the  decision  was  made  to  end  the  experiment. 
Upon  return  to  the  Monterey  Bay  Aquarium  Research 
Institute,  the  gravity  cores  were  taken  immediately  to  the 
USGS  Laboratories  in  Menlo  Park,  California  where  they 
were  refrigerated  until  the  analysis  could  be  performed  on 
May  31  -  June  1,  2000. 

Analysis  of  the  gravity  cores  was  begun  on  May  31, 
2000  at  the  USGS  Laboratories  in  Menlo  Park,  California. 
The  gravity  cores  were  sliced  into  two-centimeter  segments 
to  a  depth  of  ten  centimeters,  and  then  sliced  into  four- 
centimeter  segments.  A  Fall  Cone  Apparatus  (Model  G- 
200)  was  used  to  determine  sediment  shear  strength. 
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In  the  test,  it  is  assumed  that  the  shear  strength  of 
sediment  at  constant  penetration  of  a  cone  is  directly 
proportional  to  the  weight  of  the  cone  and  the  relation 
between  undrained  shear  strength  s  and  the  penetration  h  of 
a  cone  of  weight  Q  is  given  by: 

S(z)  =  KQ/h2,  (34) 

where  K  is  a  constant  which  depends  mainly  on  the  angle  of 
the  cone,  but  is  also  influenced  by  the  sensitivity  of  the 
clay/sediment. 


Table  1.  Physical  parameters  of  the  model  cylinders 


c 
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Mass 

(g) 

L  (cm) 

Volum 
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(cm3) 

Pm 

(g 

Ill’) 
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X 

(cm) 

Jl  (■ h ) 
(gm2) 

1 

322.5 

15.20 

191.01 

1.69 

330.5 

0.00 

( 

•087.9 

0.74 

i 

1783.0 

1.48 

( 

.233.8 

2 

254.2 

12.10 

152.05 

1.67 

271.3 

0.06 

424.6 

0.53 

' 

206.5 

1.00 

' 

312.6 

3 

215.3 

9.12 

114.61 

1.88 

235.0 

0.00 

695.2 

0.29 

577.5 

0.58 

556.8 

Table  2.  Number  of  drops  conducted  for  different 
drop  angles  and  X  -values  for  L/R  =  15/2. 


rf> 

15° 

30° 

45° 

60° 

75° 

z2 

13 

15 

15 

15 

12 

Zi 

9 

15 

15 

15 

9 

Zo 

12 

14 

15 

18 

6 

Z- 1 

0 

6 

6 

6 

0 

X-2 

2 

6 

6 

0 

0 

Four  different  cones  are  used  with  this  instrument,  each 
one  having  the  measuring  range  listed  in  Table  3.  The 
cones  are  suspended  from  a  permanent  magnet.  By 
pressing  a  knob,  the  magnet  is  moved  so  that  the  magnetic 
field  is  broken  momentarily,  and  the  cone  is  released. 
Measurements  are  taken  of  penetration  depth  and  the 
evolution  is  repeated  five  times  per  sediment  slice.  These 
values  are  then  averaged  and  correlated  with  a  table  which 
gives  shear  strength.  Previous  studies  (Chu  et  al.  2002) 
showed  that  the  sediment  parameters  are  the  most  critical 


element  in  determining  how  deep  an  object  was  buried 
when  it  came  to  rest.  During  the  experiment  at  the 
Monterey  Bay,  we  obtained  17  gravity  cores.  Sediment 
bulk  density  and  shear  strength  profiles  (Fig.  11)  show 
generally  increase  with  depth  until  approximately  6-9  cm 
below  the  water-sediment  interface. 


Table  3.  Measuring  ranges  of  the  gravity  cores 


Weight 

(g) 

Apex- 

Angle 

Penetration 

(mm) 

Undrained 
shear  strength 

(kPa) 

400 

u> 

o 

o 

4.0-15.0 

25-1.8 

100 

o 

o 

5.0-15.0 

4-0.45 

60 

60° 

5.0-15.0 

0.6  -  0.067 

10 

o 

O 

5.0-20.0 

0.10-0.0063 

Fig.  11.  Sediment  density  and  shear  strength  profiles  in 
the  Monterey  Bay  collected  during  the  cylinder  drop 
experiment  on  May  31,  2000. 

10.  Model-Data  Comparison 

The  U.S.  Navy  has  a  2D  model  in  x-z  plane 
(IMPACT28)  to  predict  cylinder’s  trajectory  and  impact 
burial.  Since  the  motion  of  cylinder  is  3D,  the  impact  burial 
prediction  using  the  2D  model  has  large  errors  (Chu  et  al., 
2000).  In  this  study,  a  new  3D  model  (called  IMPACT35) 
is  developed  on  the  base  of  momentum  balance  (7)  and 
moment  of  momentum  balance  (8)  using  triple  coordinate 
transform  ([2]  Chu  et  al.  2004)  and  cylinder  decomposition. 
To  evaluate  the  value-added  of  the  3D  model,  comparison 
among  the  observed  data  (from  Exp-1  and  Exp-2)  and 
predicted  data  using  2D  (IMPACT28)  and  3D 
(IMPACT35)  models  is  conducted. 
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10.1.  Comparison  Using  Exp-1  Data 

Improvement  from  IMPACT28  to  IMPACT35  in 
predicting  cylinders’  trajectory  and  orientation  in  the  water 
column  is  verified  using  the  Exp-1  data.  Here,  we  list  two 
examples. 

Positive  X  (Nose-Down):  Cylinder  #1  (L  =  15.20  cm, 

p  =  1.69  g  cm’3 )  with  x  =  0-74  cm  is  injected  to  the  water 
with  the  drop  angle  45°.  The  physical  parameters  of  this 
cylinder  are  given  by 
m  -  322.5  g, 

2  2  •  (35a) 

Jx  =  330.5  g  cm  ,  J2=J3  =  5783.0  g  cm 

Undersea  cameras  measure  the  initial  conditions 

xo  =  To  =  °>  zo  =  °> 

uQ  =  0,  v0  =  -1.55  m  s’1,  w0  =  -2.52  m  s’1, 

^10=°’  ^20  =  60°>  V'io  =  -95°’ 

(35b) 

col0  =  0,  co20  =  0.49  s’1,  co30  =  0.29  s’1 

Substitution  of  the  model  parameters  (35a)  and  the  initial 
conditions  (35b)  into  IMPACT28  and  IMPACT35  leads  to 
the  prediction  of  the  cylinder’s  translation  and  orientation 
that  are  compared  with  the  data  collected  during  Exp-1  at 
time  steps  (Fig.  12).  The  new  3D  model  (IMPACT35) 
simulated  trajectory  agrees  well  with  the  observed 
trajectory.  Both  show  the  same  slant- straight  pattern  and 
the  same  travel  time  (1.23  s)  for  the  cylinder  passing 
through  the  water  column.  However,  the  existing  2D 
model  (IMPACT28)  has  less  capability  to  predict  the 
cylinder’s  movement  in  the  water  column.  The  travel  time 
predicted  by  IMPACT28  is  1.5  s,  much  more  than  the 
observed  value. 


Negative  X  (Nose-Up):  Cylinder  #2  ( L  =  12.10  cm, 

p  =  1.67  g  cm’3)  with  x  =  -1-00  cm  is  injected  to  the 
water  with  the  drop  angle  30°.  The  physical  parameters  of 
this  cylinder  are  given  by 
m  -  254.2  g, 

2  2  •  (36a) 
Jx  -  271.3  gem,  J2  =  J3  =  3312.6  g  cm 

Undersea  cameras  measure  the  initial  conditions 


xo  =  To  =  °>  zo  =  °> 

uQ  -  0,  v0  =  -0.75  m  s’1,  wQ  =  -0.67  m  s’1, 


Vio =0’  y 20  ~  24°  ’  ^30  =-96°> 


®10  -  60 20  ~ 


-5.08  s 


-l 


(36b) 


a)30  =0.15  s’ 


The  predicted  cylinder’s  translation  and  orientation  are 
compared  with  the  data  collected  during  Exp- 1  at  time  steps 
(Fig.  13).  The  new  3D  model  (IMPACT35)  simulated 
trajectory  agrees  well  with  the  observed  trajectory.  Both 


show  the  same  flip-spiral  pattern  and  the  same  travel  time 
(1.73  s)  for  the  cylinder  passing  through  the  water  column. 
The  flip  occurs  at  0.1 1  s  (0.13  s)  after  cylinder  entering  the 
water  in  the  experiment  (IMPACT35).  After  the  flip,  the 
cylinder  spirals  down  to  the  bottom.  However,  the  existing 
2D  model  (IMPACT28)  does  not  predict  the  flip-spiral 
pattern. 

Eiip^riifiert  3-0  Model 


HKHjfil 


'  8 

|pn|LQ&0 

-Otfi  -ft*  -0.4  -0 1  6  0.J 


Fig.  12.  Movement  of  Cylinder  #1  ( L  =  15.20  cm, 
p  =  1.69  g  cm"3)  with  x  =  0.74  m  and  drop  angle  45° 


obtained  from  (a)  experiment,  (b)  3D  IMPACT35 
model,  and  (c)  2D  Impact28  model. 


Fig.  13.  Movement  of  Cylinder  #2  ( L  =  12.10  cm, 
p  =  1.67  g  cm  )  with  x  =  -1-00  cm  and  drop  angle  30° 
obtained  from  (a)  experiment,  (b)  3D  IMPACT35 
model,  and  (c)  2D  Impact28  model. 

10.2.  Comparison  Using  Exp-2  Data 

After  running  the  two  models  (IMPACT35  and 
IMPACT28)  for  each  gravity  core  regime  and  location,  the 
burial  depths  were  compared  with  measured  burial  depth 
data  (Fig.  14).  As  evident,  IMPACT35  improves  the 
prediction  capability.  The  existing  2D  model  (IMPACT25) 
over  predicts  actual  burial  depth  by  an  order  of  magnitude 
on  average.  However,  the  3D  model  (IMPACT35)  predicts 
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the  burial  depth  reasonably  well.  Since  the  gravity  cores 
were  taken  for  approximately  two  to  three  meters  from  the 
impact  location,  several  cores  were  taken  for  each  drop. 
This  allowed  an  average  to  be  calculated  in  order  to  yield 
more  accurate  data  for  each  drop. 


Drop  Number 

Fig.  14.  Comparison  among  observed  and  predicted 
burial  depths. 

11.  Conclusions 

(1)  A  3D  Navy’s  model  for  Mine  Impact  Burial 
Prediction  (IMPACT35)  is  developed  to  predict  the 
translation  and  orientation  of  falling  rigid  cylinder  through 
air,  water,  and  sediment.  It  contains  three  components: 
triple  coordinate  transform,  cylinder  decomposition,  and 
hydrodynamics  of  falling  rigid  object  in  a  single  medium 
(air,  water,  or  sediment)  and  in  multiple  media  (air-water 
and  water-sediment  interfaces). 

(2)  Triple  coordinate  transform  is  useful  for  modeling  the 
movement  of  rigid  body  in  air-water-sediment.  The  body 
forces  (including  buoyancy  force)  and  torques  are 
represented  in  the  E-coordinate  system,  the  hydrodynamic 
forces  (such  as  the  drag  and  lift  forces)  and  torques  are 
represented  in  the  F-coordinate,  and  the  cylinder’s  moments 
of  gyration  are  represented  in  the  M-coordinate.  The 
momentum  (moment  of  momentum)  equation  for  predicting 
the  cylinder’s  translation  velocity  (orientation)  is 
represented  in  the  E-coordinate  (M-coordinate)  system. 
Transformations  among  the  three  coordinate  systems  are 
used  to  convert  the  forcing  terms  into  E-coordinate  (M- 
coordinate)  for  the  momentum  (moment  of  momentum) 
equation. 

(3)  Cylinder  decomposition  proposed  in  this  study  leads 
to  effective  computation  of  the  forces  and  torques  as  the 
cylinder  penetrates  an  interface  between  two  media. 
Volume  decomposition  is  used  to  compute  body  forces  and 
torques.  The  surface  decomposition  is  used  to  calculate 
hydrodynamic  forces  and  torques. 


(4)  Two  cylinder  drop  experiments  were  conducted  to 
evaluate  the  3D  model.  Model-data  comparison  shows  that 
IMPACT35  improves  the  prediction  capability. 

(5)  Future  work  includes  extension  of  IMPACT35  for 
cylindrical  mines  (current  status)  to  any  shapes  of  mines, 
and  verification  of  IMPACT35  using  the  Navy’s  real-time 
mine  warfare  exercise  data  as  well  as  other  data  regarding 
the  movement  of  rigid  objects  (such  as  bombs)  in  the  air- 
water- sediment. 

Acknowledgments 

The  Office  of  Naval  Research  Marine  Geosciences 
Program  (N0001403WR20178),  Naval  Oceanographic 
Office,  and  the  Naval  Postgraduate  School  supported  this 
study. 

References 

Amone,  R.  A.,  and  Bowen,  Prediction  Model  of  the  Time  History 
Penetration  of  a  Cylinder  through  the  Air- Water- Sediment  Phases.  NCSC 
letter  report  T34,  Naval  Coastal  Systems  Center,  Panama  City,  FL,  1980. 
Boorda,  J.  M.,  "Mine  Countermeasures  -  An  Integral  Part  of  our  Strategy 
and  our  Forces."  Federation  of  American  Scientists, 
(http  ://www.fas.org/man/dod- 101  /sys/ship/weaps/docs/cnopaper.htm) 

Chu,  P.C.,  E.  Gottshall,  and  T.E.  Halwachs,  1998:  Environmental  Effects 
on  Naval  Warfare  Simulations.  Institute  of  Joint  Warfare  Analysis,  Naval 
Postgraduate  School,  Technical  Report,  NPS-IJWA-98-006,  33p. 

Chu,  P.C.,  V.I.  Taber,  and  S.D.  Haeger,  2000a:  A  Mine  Impact  Burial 
Model  Sensitivity  Study.  Institute  of  Joint  Warfare  Analysis,  Naval 
Postgraduate  School,  Technical  Report,  NPS-IJWA-00-003,  48p. 

Chu,  P.C.,  V.I.  Taber,  and  S.D.  Haeger,  2000b:  Environmental  Sensitivity 
Study  on  Mine  Impact  Burial  Prediction  Model.  Proceedings  on  the  Fourth 
International  Symposium  on  Technology  and  the  Mine  Problem,  10  pp. 
Chu,  P.C.,  T.B.  Smith,  and  S.D.  Haeger,  2001:  Mine  burial  impact 
prediction  experiment.  Institute  of  Joint  Warfare  Analysis,  Naval 
Postgraduate  School,  Technical  Report,  NPS-IJWA-0 1-007,  161p. 

Chu,  P.C.,  A.F.  Gilles,  C.W.  Fan,  J.  Lan,  and  P.  Fleischer,  2002a: 
Hydrodynamics  of  falling  cylinder  in  water  column.  Advances  in  Fluid 
Mechanics,  4,  163-181. 

Chu,  P.C.,  T.B.  Smith,  and  S.D.  Haeger,  2002b:  Mine  impact  burial 
prediction  experiment.  Journal  of  Counter-Ordnance  Technology  (CD 
Rom). 

Chu,  P.C.,  A.F.  Gilles,  C.W.  Fan,  and  P.  Fleischer,  2002c:  Hydrodynamics 
of  falling  mine  in  water  column.  Journal  of  Counter-Ordnance  Technology 
(CD  Rom). 

Chu,  P.C.,  C.W.  Fan,  A.D.  Evans,  and  A.  Gilles,  2003:  Triple  coordinate 
transforms  for  prediction  of  falling  cylinder  through  the  water  column. 
Journal  of  Applied  Mechanics,  in  press. 

Hurst,  R.B.  Mine  Impact  Burial  Prediction  Model  -  Technical  Description 
of  Recent  Changes  and  Developments.  Defense  Scientific  Establishment, 
Auckland,  New  Zealand,  Report  149. 

Lehr,  S.E.,  Mine  Warfare:  To  Enable  Maneuver,  Amphibious  Warfare 
Conference,  26  April  2000.  (http:// 

www.exwar.org/awcfinal/6th/n852mines.ppt) 

Lott,  D.F.,  K.  Williams,  and  D.  Jackson,  Mine  Burial  in  Carbonate 
Sediments.  Proc.  of  the  Technology  and  Mine  Problem  Symposium, 
November  1996,  Naval  Postgraduate  School. 

Rhodes,  J.E.,  G.  Holder,  Concept  for  Future  Naval  Mine  Countermeasures 
in  Littoral  Power  Projection,  1  May  1998.  (http:// 
192. 156.75. 102/mcm.htm) 

Satkowiak,  L.  J.,  User's  Guide  for  the  Modified  Impact  Burial  Prediction 
Model.  NCSC  TN  884-87.  Naval  Coastal  Systems  Center,  Panama  City,  FI, 
1987. 


10 


